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Taking into account of entropy generation during evolution of a viscous fluid, we have estimated 
inverse Knudsen number, ideal hydrodynamic limit for elliptic flow and QGP viscosity to entropy 
ratio in y^sjvjv=62 and 200 GeV Cu+Cu/Au+Au collisions. Viscosity to entropy ratio is estimated 
as rj/s — 0.17 ± 0.10 ± 0.20, the flrst error is statistical, the second one is systematic. In a central 
Au+Au collision, inverse Knudsen number is ~ 2.80 ± 1.63, which presumably small for complete 
equilibration. In peripheral collisions it is even less. Ideal hydrodynamic limit for elliptic flow is 
~40% more than the experimental flow in a central collision. 



o 

(N 

^■ 
CD 

(N 



> 
00 

o 
o 



PACS numbers: 47.75.+f, 25.75.-q, 25.75.Ld 



Recent experiments at RHIC produced convincing ev- 
idences that a collective matter is created in Au-|-Au col- 
lisions [H-Q- The evidences come mainly from observing 
finite elliptic flow in non-central collisions, which is now 
regarded as a definitive signature of collective effect [IjQ. 
Elliptic flow measure the azimuthal correlation of pro- 
duced particles with respect to the reaction plane. It is 
also best understood in a collective model like hydrody- 
namics i?! . In a non-central collision, the reaction zone is 
spatially asymmetric. Differential pressure gradient con- 
vert the spatial asymmetry in to momentum asymmetry. 
Ideal hydrodynamics has been quite successful explaining 
a large part of the experimental elliptic flow How- 
ever it is now realized that the experimentally measured 
scaling of integrated V2 with multiplicity or with collision 
centrality is not in agreement with ideal hydrodynamics. 
While ideal hydrodynamics predicts approximate scaling 
[8|, in ex per iments scaling is violated 

UM- As dis- 
cussed in [III , violation of the scaling can be understood 
as an indication of incomplete thermalization. Ideal hy- 
drodynamics require local thermal equilibration. Devia- 
tion from the local equilibrium can lead to a characteris- 
tic dependence of the eccentricity scaled elliptic flow on 
charged particles multiplicity. Degree of thermalization 
in the fluid produced in Au-|-Au collisions can be charac- 
terized by the dimensionless parameter, Knudsen number 
(K) 11]. By definition, inverse of the Knudsen number 
is the number of collisions per particles. 
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where R is the characteristic size of the system, n is the 
particle density and cr is the inter particle cross section. 



Validity of hydrodynamics require that K 



» 1, so 



that large number of collisions can bring the system to 
local thermal equilibrium. The opposite limit, K^^ << 



1 is the Knudsen regime where hydrodynamics become 
inapplicable. The simple formula, 
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proposed in [ll| give qualitatively correct behavior of the 
experimental elliptic flow. In the limit of small Knudsen 
number, experimental flow approaches the ideal hydro- 
dynamic limit with a correction linear m. K. In 
the other extreme limit of large flow is proportional 
to the Knudsen number K. In EqlJl K^^ is a number 
of the order of unity, whose precise value can be deter- 
mined only from explicit transport calculations. From 
Monte-Carlo simulation of transport equations ATo was 
estimated, A'o = 0.70 ± 0.03 [12]. In [11] it was argued 
that inverse of the Knudsen number can be deter- 
mined from the experimental data, as it is proportional 



to 



1 dN 



~S~dy^ where ^ is the total multiplicity density and S 
is a measure of the transverse area of the collision zone, 
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In Eq[21 Cs is the speed of sound of the medium and a 
is the intcr-particle cross section. Eql2]and[3]connect two 
experimental observables, elliptic flow and particle multi- 
plicity and can be used to determine unknown quantities 
e.g. ideal hydrodynamic limit of elliptic fiow (■^)' , the 
combination of parameters KquCs- 

However, there is a serious flaw in the derivation of 
Eq|3] and the conclusions derived from it could be mis- 
leading. Eq|3]was obtained with the assumption that the 
total particle number is conserved throughout the evolu- 
tion 11]. The assumption is justified in an isentropic 
expansion, i.e. one dimensional evolution of ideal fiuid, 
when entropy density (s) times the proper time (r) is a 



constant. Under such condition, 
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However, in a viscous evolution, entropy is generated and 
initial and final state entropy are not same and the as- 
sumption is clearly violated. Only in systems with very 
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FIG. 1: (color online) PHOBOS data for the centrality de- 
pendence of eccentricity scaled elliptic flow in Cu+Cu and 
Au+Au collisions at ^/sNN— 62 and 200 GeV. The solid line 

is the fit to the data with (i^)'''=0.33 and (^:p^^^^0.83. 
The dashed line is the fit in the ideal fluid approximation, 
(f)"'= 0.19 and (-^f)=0. 



small viscosity, the assumption may be approximately 
valid, but not in systems where sufficient entropy is gen- 
erated. Explicit numerical simulations indicate that in 
Au+Au collisions, entropy generation can be substan- 
tial, e.g. ~ 20%, 30% and 50% in fluid evolution with 
viscosity to entropy ratio, 77/3=0. 08, 0.12 and 0.16 [l3 |. 
One may argue that unlike in explicit numerical simula- 
tion of viscous hydrodynamics, Knudsen number ansatz 
does not require the entire dynamic range of evolution. 
Time of validity of Knudsen number ansatz is t < R/cs 
[Tl| . For Cs = -y/l/S, and characteristic size R « 1-2 fm, 
the time scale is r <1.7-3.5 fm. Though the time scale 
is small compared to the entire dynamic range of evolu- 
tion, it is large enough for significant entropy generation. 
In explicit simulation of viscous hydrodynamics, entropy 
generation is fast and most of the entropy is generated 
with first 2-4 fm of evolution [l^ . Thus even though 
Knudsen number ansatz does not require full dynamic 
range of evolution, in the time scale for validity of Knud- 
sen ansatz most of the entropy will be generated. 

In the present paper accounting for the entropy gen- 
eration in viscous evolution, we generalise Eql31 As in 
PH . we also consider one dimensional Bjorken longitudi- 
nal expansion. If rj and C is the shear and bulk viscosity 
coefficients, for Bjorken flow, energy-momentum conser- 
vation equation, dfj.T^" = reduces to [TgI. [T7|. 
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where s is the entropy density and r is the proper time. 
In general bulk viscosity is much less than the shear vis- 
cosity, however, near the phase transition region bulk 
viscosity can be large [3 In the following we ne- 
glect the bulk viscosity. We further assume that shear 
viscosity is proportional to cube of the temperature. For 



r] ocT^, Eq|4]can be analytically integrated between time 
Ti and Tf. For TV >> r^, the final state entropy can be 
written as [l6l.[l7l|. 
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In EqlSl Ti is the temperature at the time scale 
Equating final state entropy density s f with particle mul- 
tiplicity per unit transverse area {S), TfSf w 3.6^^ 
[isj . we obtain the particle density at the time scale Tj 



as. 
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Inserting Eql6] in EqlU in the time scale ti < R/cs, 
i.e. as long as the transverse size of the system does not 
vary significantly [Tl| , inverse of Knudsen number can be 
obtained as. 
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One immediately observes that neglect of entropy gen- 
eration during evolution will over estimate K~^, by the 



factor 



1 
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As it will be shown below, ex- 
perimental data indicate that the factor could be large, 
~ 2-7. 

Inserting Eql7]in Eql^l eccentricity scaled elliptic flow 
now can be related to observed particle multiplicity as. 
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Experimental data on elliptic flow and particle mul- 
tiplicity can be fitted with Eql8] to obtain estimates of 

the hydrodynamic limit of elliptic flow (^)' and vis- 
cosity to entropy ratio in unit of initial time and tem- 
perature (^f^ EqlHalso involve the quantity KoaCs- 

Kq, a and Cg are known approximately. For example, 
inter parton cross section is expected to be small, tT=3-4 
mb. From transport calculations, Kq was estimated as 
Kq — 0.7 ± 0.3. The speed of sound of QGP medium is 
expected to be Cg ~ In the following, we flx a=3 

mb, and Cg = \/l/3. 

In FigE PHOBOS ^- 22.] data for the centrality de- 
pendence of (participant) eccentricity scaled charged par- 
ticles elliptic flow, in Au+Au and Cu+Cu collisions at 
77)^=62 and 200 GeV are shown. PHOBOS data 
have large error bars and within the error data do not 
show any system size or energy dependence. All the 
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FIG. 2: Viscosity to entropy ratio r;/s as a function of the 
parameter from the analysis in [l^. The sohd hne 

is a straight hne fit. 



four data sets are fitted together with Eq|51 The solid 



Data are well explained, 

, ih 



line in Fig IT] shows the fit. 

value is also small, x^/N w 0.1. However, (2^)"' and 
{^■:pjr^ can be determined only with large uncertain- 
ties, (2^)''' = 0.33±0.12, i -^A = 0.83±0.51. Uncer- 



V 



would be reduced with bet- 



tainty in (^22.^*'' and 
ter quality data. While definitive conclusions cannot be 

made due to large uncertainty in (^) * , the central value 
suggests that even in central Au+Au collisions, ideal hy- 
drodynamic limit is not reached. In central/mid-central 
collisions, experimental flow is only ^ 60% of the ideal 
fluid limit. In peripheral collisions, it is even less. For 
comparison purpose, in Fig HI we have shown the flt to 

the data in the ideal fluid approximation, ^7:^7^ = 0. 
In the ideal fluid approximation ideal hydrodynamic limit 
for elliptic flow is estimated as, (^a)*'' = 0.19 ± 0.04. 
However, description to the data is comparatively poor, 
X^/ ~ 1.3 is ~ 10 times larger than that obtained in 
viscous evolution. 

With the estimate of (yyyjr^ , we can compute inverse 
Knudsen number from EqjT] In tablelH we have listed 
for different centrality ranges of Au+Au collisions. 
Particle multiplicities are taken from the PHOBOS data 
[20l-l22'| . For comparison, we have also listed the values in 
the ideal fluid approximation. decreases by a factor 
of 2-7 when entropy generation is accounted for. In the 
ideal fluid approximation, « 8 in a central collision. 
It is reduced to ~2.8 if entropy generation during evo- 
lution is accounted for. While it is debatable whether 
ss 8 can lead to complete equilibration, it is un- 
likely that complete thermalization will be achieved with 
« 2.8. In more peripheral collisions, equilibration 
is certainly incomplete. 

QGP viscosity is an important parameter charac- 
terizing QGP medium. String theory based models 
(ADS/CFT) give a lower bound on viscosity to entropy 



TABLE I: Inverse Knudsen number {K~^), in ideal and vis- 
cous fluid, as a function of collision centrality in Au-|-Au col- 
lision. Also listed are charged particles multiplicity per unit 
transverse area ( 



1 dJVcfc - 

S dy 



from PHOBOS experiment [20|-l2 
and characteristic size R of the system. Inter parton cro ss se c- 
tion is assumed to be cr= 3 mb, speed of sound Cs = \/T/3. 
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(AjtA =0.83 ±0.51 




centrality(%) 




(fm) 






0-3 


30.46 


2.01 


2.80 ± 1.63 


7.91 


3-6 


28.84 


1.93 


2.65 ±1.55 


7.49 


6-10 


27.06 


1.84 


2.49 ± 1.45 


7.03 


10-15 


24.95 


1.75 


2.29 ± 1.33 


6.48 


15-20 


23.17 


1.64 


2.13 ± 1.24 


6.02 


20-25 


21.39 


1.55 


1.97 ± 1.15 


5.56 


25-30 


19.44 


1.45 


1.79 ± 1.04 


5.05 


30-35 


17.50 


1.38 


1.61 ±0.94 


4.55 


35-40 


14.26 


1.30 


1.31 ±0.76 


3.70 


40-45 


11.83 


1.23 


1.09 ±0.63 


3.07 


45-50 


10.69 


1.18 


0.98 ±0.57 


2.78 



TABLE II: Listed are some estimates of QGP viscosity to 
entropy ratio from experimental data in Au+Au collisions at 
RHIC. The observables analyzed are also listed. 



SI. no. 




Experimental observable 


1 


0.88 ± 0.38 ± 1.76 


<j) meson's (A''), (pr) and V2 fl4] 


2 


1.0-3.77 


Pt fluctuations [26] 


3 


1.4-2.4 


V2 scaling violation [27] 


4 


«2 


V2 scaling violation [2^ 


5 


1.13 ± 0.19 ± 1.26 


V2 scaling violation [2£] 


6 


1.51 ±0.38 


V2 scaling violation [30] 


7 


1.3-2.0 


heavy quark energy loss [3L] 


8 


1.45 ±0.06 


V2 [32] 


9 


< 1.51 


Pt spectra of tt, A" and p [33] 


10 


2.14 ± 1.26 ±2.56" 


"^'^ S dy 



"present work 



ratio of any matter, inrj/s > 1 [23|. In [2J, |25|, from ex- 
perimental data, a phenomenological upper bound was 
conjectured, Awrj/s < 5. We have obtained viscosity to 
entropy ratio in unit of initial time and temperature. It 
can be converted to more comprehensible viscosity to en- 
tropy ratio if the initial time and temperature scale is 
known. Recently, in [T3] STAR data on (f> mesons multi- 
plicity, mean pT and integrated V2 were analyzed in ideal 
and viscous fluid dynamics. At the initial time Ti— 0.6 
fm, ideal or viscous fluid was initialized to reproduce ex- 
perimental (f> meson multiplicity. Viscous fluid requires 
less initial temperature than an ideal fluid. Results of 
the analysis are shown in Figj^l where rj/s as a func- 
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tion of the parameter 




is shown. The sohd hne 



in FigHis a straight hne fit, 'q/s = 0.2 (7^7) • The 

relation can be used to convert extracted (^7^7^ to 

'q/s. We obtain, -q/s = 0.17 ± 0.10. From the PHO- 
BOS data, viscosity to entropy ratio can be determined 
only within ~ 60% accuracy. Evidently, much better 
quality data are required for more precise determination 
of viscosity to entropy ratio. As mentioned earlier, the 
estimate was obtained with Kq = 0.7, = and 
a = 0.3/m^ corresponding to, Kqo-Cs ~ 0.121. Estimate 
of viscosity depend on the value of KoaCs- Systematic 
uncertainty in -q/s is ~ 120% due to a factor of 2 un- 
certainty in KqctCs- We then estimate QGP viscosity to 
entropy ratio as -q/s = 0.17 ± 0.10 ± 0.20, the first error 
is statistical, the second is systematic. Systematic un- 
certainty will increase if uncertainty in initial time and 
temperature scale is included. The estimated value is 
weh with the two bounds, 1 < An-q/s < 5 [IS-il. In 
tableini present estimate for QGP viscoisty is compared 
with some recent estimates. One may note that presently 
estimated rj/s is similar to the values obtained in previous 
extractions (27|,^i, which disregarded entropy genera- 
tion in the Knudsen ansatz. The reason is understood. 
Experimental data include the effect of entropy gener- 
ation. Thus even if viscous effects are neglected in the 
Knudsen ansatz, the fitted paratmeters ' and Koacg 



(see EqlH]) will include the effect. 

We note that the present estimate should be consid- 
ered as an upper limit for QGP viscosity. We have ne- 
glected bulk viscosity. Experimental data include the 
effect of bulk viscosity. Neglect of bulk viscosity will 
be compensated by increasing rj/s. Also, we have ne- 
glected transverse expansion. Experimental data also in- 
clude the effect of transverse expansion. One observes 
from Eql7]that will decrease if transverse expansion 
is included (transverse area of system at freeze-out will 
be larger than the initial area). Neglect of the transverse 
expansion will be compensated again by increasing rj/s. 

To conclude, taking into account that entropy is gen- 
erated during evolution of a viscous fluid, we have gener- 
alized a relation between inverse Knudsen number 
and particle density PHOBOS data on the cen- 

trality dependence of elliptic flow indicate that K^^ is 
overestimated by a factor ^ 2-7 if entropy generation is 
neglected. We have also estimated ideal hydrodynamic 
limit for elliptic flow and QGP viscosity to en- 

tropy ratio 77/s. Estimated (^^y'^ is ^40% larger than 
the scaled flow in a central collision. Estimated viscos- 
ity to entropy ratio, Anrj/s = 2.14 ± 1.26 is well within 
the ADS/CFT lower bound and phenomenological upper 
bound, 1 < Anrj/s < 5 . 



BRAHMS Collaboration, I. Arsene et al, Nucl. Phys. A 

757, 1 (2005). [17] 

PHOBOS Collaboration, B. B. Back et al, Nucl. Phys. 

A 757, 28 (2005). [18 

PHENIX Collaboration, K. Adcox et al., Nucl. Phys. A [19 

757 184 (2005). 

STAR Collaboration, J. Adams et ai, Nucl. Phys. A 757 [20 
102 (2005). 

J. Y. OUitrault, Phys. Rev. D 46, 229 (1992). [21 
A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58, 
1671 (1998) [22 
P. F. Kolb and U. Heinz, in Quark-Gluon Plasma 3, [23 
edited by R. C. Hwa and X.-N. Wang (World Scientific, 
Singapore, 2004), p. 634. [24 
H. Song and U. W. Heinz, Phys. Rev. C 78, 024902 
(2008). [25 
S. A. Voloshin [STAR Collaboration] , J. Phys. G 34, S883 [26 
(2007) . 

S. A. Voloshin [STAR Collaboration], AIP Conf. Proc. [27 
870, 691 (2006). 
R. S. Bhalerao, J. P. Blaizot, N. Borghini and J. Y. 01- [28 
litrault, Phys. Lett. B 627, 49 (2005). 

C. Gombeaud and J. Y. OUitrault, Phys. Rev. C 77, [29 
054904 (2008) [30 
R. C. Hwa and K. Kajantie, Phys. Rev. D 32, 1109 [31 
(1985). 

A. K. Chaudhuri, Phys. Lett. B 681, 418 (2009). [32 
A. K. Chaudhuri, Phys. Rev. C 74, 044904 (2006) [33 
P. Danielewicz and M. Gyulassy, Phys. Rev. D 31, 53 



(1985). 

A. Muronga, Phys. Rev. Lett. 88, 062302 (2002) 
[Erratum-ibid. 89, 159901 (2002)] 

D. Kharzeev and K. Tuchin, JHEP 0809, 093 (2008) 

F. Karsch, D. Kharzeev and K. Tuchin, Phys. Lett. B 
663, 217 (2008) 

B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. 
C 72, 051901 (2005) 

B. Alver et al. [PHOBOS Collaboration], Phys. Rev. 
Lett. 98, 242302 (2007) 

B. Alver et al, Phys. Rev. C 77, 014906 (2008) 

G. Policastro, D. T. Son and A. O. Starinets, Phys. Rev. 
Lett. 87, 081601 (2001). 

M. Luzum and P. Romatschke, Phys. Rev. C 78, 034915 
(2008). 

H. Song and U. W. Heinz, J. Phys. G 36, 064033 (2009). 
S. Gavin and M. Abdel-Aziz, Phys. Rev. Lett. 97, 162302 
(2006). 

H. J. Drescher, A. Dumitru, C. Gombeaud and J. Y. 01- 
litrault, Phys. Rev. C 76, 024905 (2007). 
H. Masui, J. Y. OUitrault, R. Snellings and A. Tang, 
Nucl. Phys. A 830, 463C (2009). 

R. A. Lacey et al, Phys. Rev. Lett. 98, 092301 (2007). 
A. K. Chaudhuri, Phys. Rev. C 81, 044905 (2010). 
A. Adare et al. [PHENIX CoUaboration], Phys. Rev. 
Lett. 98, 172301 (2007). 

A. K. Chaudhuri, J. Phys. G 37, 075011 (2010). 

V. Roy and A. K. Chaudhuri, a rXiv:1003.119 5 [nucl-th]. 



